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Abstract-  To  date  various  image  filtering  approaches  are  re¬ 
ported,  most  of  them  enhance  the  signal-to-noise  ratio(SNR) 
in  different  levels  with  loss  of  some  useful  information.  In 
this  paper,  an  adaptive  template  filtering  method  is  described 
which  can  be  used  to  increase  SNR  and  keep  the  important 
edge  information  of  medical  images.  In  the  proposed  method, 
for  each  pixel,  an  optimal  template  of  different  shape  and  co¬ 
efficients  is  selected  automatically  depending  on  its  neighbor 
pixels.  Not  being  simple  extension  of  the  2D  method,  se¬ 
quential  filtering  of  slices  in  three  orthogonal  directions  imple¬ 
mented  the  3D  algorithm,  which  improved  existing  adaptive 
template  filtering  technique,  avoiding  a  burden  of  enormous 
search.  Simulation  and  MRI  image  tests  ,both  2D  and  3D, 
show  that  the  new  adaptive  template  filter  provides  higher 
SNR  and  sharper  edges. 

Keywords:  Adaptive  template  filtering,  magnetic  reso¬ 
nance  imaging,  edge  extraction 

1.  INTRODUCTION 

In  medical  imaging,  if  the  SNR  is  too  small  or  the  con¬ 
trast  too  low  it  becomes  very  difficult  to  detect  anatomical 
structures  because  tissue  characterization  fails.  For  a  vi¬ 
sual  analysis  of  medical  images,  the  clarity  of  details  and 
the  object  visibility  are  important,  whereas  for  image  pro¬ 
cessing  a  high  SNR  is  required  because  most  of  the  image 
segmentation  algorithms  are  very  sensitive  to  noise.  A 
major  obstacle  to  segmentation  by  multiple  thresholding 
or  by  multivariate  statistical  classification  is  insufficient 
data  quality. 

When  developing  a  filtering  method  for  medical  image 
data,  image  degradation  by  blurring  or  by  artifacts  re¬ 
sulting  from  a  filtering  scheme  is  not  acceptable.  The 
following  requirements  should  ideally  be  fulfilled1: 

a.  minimize  information  loss  by  preserving  object  bound¬ 
aries  and  detailed  structures. 

b.  Efficiently  remove  noise  in  region  of  homogeneous 
physical  properties,  and 

c.  Enhance  morphological  definition  by  sharpening  dis¬ 
continuities. 

To  date  various  image  filtering  approaches  are  reported, 
most  of  them  enhance  the  signal-to-noise  ratio(SNR)  in 
different  levels  with  loss  of  some  useful  information.  In 
this  paper,  an  adaptive  template  filtering  method  is  de- 
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important  edge  information  of  medical  images.  In  the 
proposed  method,  for  each  pixel,  an  optimal  template  of 
different  shape  and  coefficients  is  selected  automatically 
depending  on  its  neighbor  pixels. 

2.  METHOD 

In  linear  spatial  filtering,  the  content  of  a  pixel  is  given  the 
value  of  the  average  brightness  of  its  immediate  neighbors. 
Simple  spatial  averaging,  often  called  ”  low-pass  filtering” , 
does  reduce  the  amplitude  of  noise  fluctuations,  but  also 
degrades  sharp  details  such  as  lines  or  edges.  The  filtering 
does  not  respect  region  boundaries  or  small  structures, 
and  the  resulting  images  appear  blurry  and  diffused.  This 
undesirable  effect  can  be  reduced  or  avoided  by  the  de¬ 
sign  of  nonlinear  filters2.  In  digital  image,  edges  are  often 
relative  with  the  local  intensity  gradient  of  the  images, 
and  can  be  extracted  with  the  intensity  gradient  meth¬ 
ods.  The  operator  of  intensity  gradient  mainly  retains 
horizontal  and  vertical  edges,  whereas  Roberts’  operator 
can  detect  slanting  edges  in  45  or  135  degrees  with  the 
horizontal  direction. 

Another  approach  is  adaptive  filtering,  which  entails 
a  tradeoff  between  smoothing  efficiency,  preservation  of 
discontinuities,  and  the  generation  of  artifacts. 

Generic  images  can  be  supposed  as  composing  of  some  2D 
or  3D  limited  regions,  which  are  homogeneous  or  slow- 
changeable,  and  seperated  by  discontinuous  boundaries. 
The  core  of  adaptive  template  filtering  is  to  search  a  ap¬ 
propriate  template,  which  best  matches  the  limited  ho¬ 
mogeneous  region  enclosing  the  considering  pixel  .  Using 
the  template, the  best  filtering  results  can  be  obtained. 
The  algorithm  defined  a  set  of  templates.  For  a  given 
pixel,  its  neighbor  pixel  may  be  classified  as  object  or 
background.  Totally,  we  have  En=2Ug  1  =  255  tem¬ 
plates.  For  the  selection  of  the  best  template,  the  follow¬ 
ing  two  values  must  be  calculated: 


where  Nj  is  the  number  of  non-background  pixels  in 
template  j ,  rrij  is  the  intensity  mean  of  pixels  in  template 
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j,  and  <jj  is  the  standard  deviation. 

For  MRI  images,  intensity  variations  in  the  limited 
homogeneous  region  of  a  template  are  mainly  caused  by 
random  noise.  Since  noise  of  MRI  image  is  independent 
of  signal,  the  standard  deviation  is  relatively  small  in  flat 
template.  For  this  region,  a  template  should  be  chosen 
which  has  most  number  of  non-background  neighbor  pix¬ 
els,  so  that  maximum  noise  suppression  can  be  obtained. 


3.  RESULTS 

To  evaluate  the  algorithm  we  generate  a  simulated  image 
on  which  there  are  some  concentric  circles  and  radiation 
lines(Fig.l,  top  left).  With  this  image  we  can  evaluate  the 
algorithm  by  its  ability  of  keeping  edges  in  different  direc¬ 
tions.  First  we  add  Gaussian  noise  on  the  image(Fig.l, 
bottom  left).  Then  two  conventional  filtering  methods, 
Intensity  gradient  method  and  Roberts’  filter,  and  the 
adaptive  template  filtering  are  applied  to  the  image. 


T  =  arg max Nj {Tj \aj  <t},  (3) 

1 3 

Whereas  in  the  edge  region,  the  sharper  change  of  in¬ 
tensity  makes  the  standard  deviation  larger.  An  edge 
template  with  minimum  standard  deviation  should  be 
used: 

T  =  arg  min  aj{Tj  \cTj  >  t},  (4) 

1 3 

The  use  of  a  flat  template  or  an  edge  template  is  de¬ 
termined  by  a  threshold  value  t.  In  order  to  speed  up 
the  selection  of  the  best  template, we  may  start  the  calcu¬ 
lation  of  the  standard  deviation  for  the  template  having 
most  non-background  pixels,  then  templates  with  lesser 
non- background  pixels.  The  first  found  template  which 
standard  deviation  is  less  than  threshold  value  t  is  the 
best  template.  If  all  template’s  standard  deviation  are 
greater  than  t,  then  the  template  with  minimum  stan¬ 
dard  deviation  is  used. 

Choosing  the  best  template, the  filter  coefficients  can  be 
obtained  by  2D  adaptive  least  square  error3: 


T.r 


a2k{x,y)I{x,y)  +  crlmj 

(TUx,y)  +  (Tl 


(5) 


°k(x’y)  =  nmx(0 ,a2(x,y)  -  a2),  (6) 

where  Tx_y  is  the  filtered  output.  a2(x,y)  and  of,  are 
input  image  variance  at  point  (x,y)  and  noise  variance  of 
the  image  respectively. 

Principally  above  mentioned  method  of  2D  adaptive 
filtering  can  be  extended  to  3D  volume  data  application. 
But  in  3D  case,  the  total  number  of  possible  templates 
is  X]nL2  ^V1  •  The  exhaustive  search  of  the  templates 
is  almost  impossible  with  current  computers.  We  do  not 
use  the  3D  templates  for  the  volume  dataset.  Actually, 
for  a  given  voxel,  the  2D  adaptive  template  filtering  are 
performed  in  three  orthogonal  planes  through  the  voxel 
sequentially.  For  MRI  volume  data,  the  whole  brain  can 
be  thought  of  as  a  stacked  slices.  The  three  orthogo¬ 
nal  planes  are  axial,  cortical  and  sagittal  slices.  First, 
the  slices  in  axial  planes  are  filtered  with  the  2D  adap¬ 
tive  template  filtering  algorithm,  then  filtered  in  cortical 
planes, finally  in  segit.tal  planes.  We  get  good  results  of 
3D  filtering  avoiding  enormous  computations. 


filtered  results 


Figure  1.  Edge  extraction  of  adaptive  filtering  and  con¬ 
ventional  methods. 
top:extracted  lines, 
bottonusubtracted  images. 


The  comparison  of  the  methods  are  shown  in  the  top 
right  of  fig.  1 ,  from  which  we  can  see  that  the  two  con¬ 
ventional  methods  can  only  retain  edges  in  particular  di¬ 
rections,  but  adaptive  template  filtering  keep  the  edges 
well  in  all  directions.  The  lower  right  of  fig.  1  shows  the 
difference  images  of  the  filtered  and  original  images.  Ob¬ 
viously,  the  proposed  method  has  best  filtering  effect. 

For  the  validation  of  adaptive  template  filtering  al¬ 
gorithm  with  medical  images,  the  MRI  image  data  were 
downloaded  from  MNI’s  BrainWeb4,  which  has  181  slices, 
slice  thick  is  1mm.  One  healthy  subject  was  scanned  27 
times,  and  the  acquired  data  were  summed  and  averaged. 
The  resulting  image  has  very  high  SNR(fig.2,(a)).  To 
test,  the  performance  of  the  filter,  some  Gaussian  noise 
was  added  (fig.2,(b)).  Fig.2,(c)  shows  that  after  adaptive 
template  filtering,  even  the  edges  of  some  small  anatom¬ 
ical  structures  can  be  seen  clearly  on  the  resulting  im¬ 
age.  The  proper  a  value  is  also  in  the  range  of  0.8-1. 2. 
Fig.2,(d)  shows  the  difference  image  of  images  after  and 
before  filtering. 
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Figure  2.  Filtered  results  of  MRI  image,  a. Original  im¬ 
age,  b. Noisy  image,  c. After  adaptive  filtering,  d.  Sub¬ 
tracted  image  X  10. 


Similar  as  in  2D  test,  the  simulated  3D  volume  data 
constructs  a  sphere  with  thin  surface.  Originated  from 
the  sphere  center  there  are  many  symmetrical  radiation 
lines.  Before  filtering  some  random  noise  was  added. 
Then  the  volume  data  were  filtered  slice  by  slice  in  three 
orthogonal  directions  sequentially  using  2D  adaptive  tem¬ 
plate  filtering  algorithm. 

We  have  also  done  the  experiments  for  MNI  3D  vol¬ 
ume  data  using  this  3D  adaptive  template  filtering  algo¬ 
rithm,  and  obtained  good  filtering  results. 

4.  CONCLUSION 

A  novel  approach  to  the  noise  suppression  and  edge  ex¬ 
traction  in  medical  images  has  been  developed  and  tested. 
Unlike  conventional  filtering,  where  the  template  shape 
and  coefficients  are  fixed,  multiple  templates  are  defined 
in  the  proposed  algorithm.  For  each  pixel,  an  optimal 
template  is  selected  automatically  depending  on  its  neigh¬ 
bor  pixels. The  adaptive  template  implements  best  filter¬ 
ing  of  medical  images.  Simulation  and  MRI  image  tests 
,both  2D  and  3D,  show  that  with  the  new  adaptive  tem¬ 
plate  filter  the  important  edge  information  of  medical  im¬ 
ages  can  be  retained  while  eliminating  noise.  Not  being 
simple  extension  of  the  2D  method,  sequential  filtering 
of  slices  in  three  orthogonal  directions  improved  existing 
adaptive  template  filtering  technique,  avoiding  a  burden 
of  enormous  search. 


